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Abstract. We show that harmonic vibrations in amorphous silicon can be decomposed to 
transverse and longitudinal components in all frequency range even in the absence of the well 
defined wave vector q. For this purpose we define the transverse component of the eigenvector 
with given oj as a component, which does not change the volumes of Voronoi cells around 
atoms. The longitudinal component is the remaining orthogonal component. We have found 
the longitudinal and transverse components of the vibrational density of states for numerical 
model of amorphous silicon. The vibrations are mostly transverse below 7 THz and above 15 
THz. In the frequency interval in between the vibrations have a longitudinal nature. Just this 
sudden transformation of vibrations at 7 THz from almost transverse to almost longitudinal 
ones explains the prominent peak in the diffusivity of the amorphous silicon just above 7 THz. 


1. Introduction 

The propagation of vibrational excitations in disordered media is one of the important problems 
in condensed matter physics. Microscopic nature of these vibrational excitations is still poorly 
understood, despite that they are responsible for the specific heat, the thermal conductivity and 
the sound propagation in amorphous materials. 

In 1999 Allen and Feldman proposed a simple classification of vibrations in disordered 
media [T] , which was based on the numerical simulations of amorphous silicon. The low-frequency 
vibrations are longitudinal (LA) or transverse (TA) plane waves (acoustic phonons). The mean 
free path of phonons decreases rapidly with frequency due to structural or dynamical disorder. 
At some frequency cjir ~ 4 THz the mean free path of phonons becomes comparable with its 
wavelength. It is a so-called loffe-Regel criterion. However, the mobility edge wioc for vibration 
localization is much higher then wir what is not the case for electron localization. Vibrations 
in amorphous materials in this relatively wide frequency range wir, < uj < tuioc are delocalized 
but not propagated as phonons. These vibrations were called diffusons because they spread by 
means of diffusion mm- 

However, this picture cannot clearly explain the prominent peak (more than 3 times in the 
value) in the diffusivity of vibrations in amorphous silicon between 7 and 13 THz [U [3]. Allen 
and Feldman noticed that the sharp rise of the diffusivity at 7 THz corresponds to the end 
of TA branch in crystalline Si which is responsible for the local minimum in the vibrational 
density of states. So the frequency 7 THz marks a point in the spectrum where vibrations 
change their character from somewhat TA-like, to somewhat LA-like with much larger group 


velocity. However the notions of transverse and longitudinal vibrations are ill defined for such 
big frequencies in amorphous silicon. 

In this paper we will generalize the notions of transverse and longitudinal vibrations for 
disordered systems like amorphous silicon. As an example we study numerically the model 
amorphous silicon (a-Si) system consisting of = 32768 atoms contained in a periodic cubic 
box of lengths Lx = Ly = of approximately 87 A. The technical details of the preparation of 
the a-Si sample have already been presented in Ref. [3]. The Si-Si interaction in the amorphous 
silicon studied here is well described by the Stillinger-Weber potential [5]. 


2. Longitudinal and transverse components of displacements 

Low-frequency vibrations even in amorphous media are well-defined plane longitudinal and 
transverse waves. These vibrations can be described in continuous medium approximation. 
The displacement fields u(r) for longitudinal (L) and transverse (T) waves have a form 


UL,T(r) = u[°^exp(fqr), 
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However, this definition contains the phonon wavevector q, which is ill-defined for high-frequency 
vibrations in amorphous medium. To generalize this approach let us note that 


divu^ = iqu.j, = 0. 


( 3 ) 


The transverse displacement UT(r) has zero divergence so it conserves local volumes. 

The natural analog of the local volumes in amorphous media are Voronoi cells around each 
atom. The Voronoi cell Vi, associated with the atom i is the set of all points in the space around 
this atom whose distance to the atom position r* is not greater than their distance to the other 
atoms Tj [6]. This type of cells is also known as Wigner-Seitz cells in crystallography for regular 
crystals. Figure shows a simple example of Voronoi cells for longitudinal and transverse waves 
in a simple quadratic lattice. A longitudinal wave evidently changes the volumes of Voronoi 
cells while the transverse wave does not change the volumes despite the change in the form of 
the cells. 

A displacement of atoms in amorphous media may (or may not) change volumes of Voronoi 
cells. We will call the displacement of atoms Uj to be transverse if it does not change the volumes 
of all Voronoi cells. For that let us introduce a matrix A which is responsible for the relative 
change of the fth Voronoi cell volume V) under jth atom displacement in the direction a 
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The explicit formula for the matrix A will be derived in the next section. Using this matrix the 
displacement of jth atom in the direction a, uja produces the following relative change of the 
Voronoi cell volumes 




( 5 ) 


ja 


In the matrix notation this equation reads £ = Au where A is a rectangular N x 3V matrix 
(with N being the number of particles) and tt is a displacement vector with 2>N elements. The 
matrix A is a discrete analog of the divergence operator. 

By definition the transverse component of the displacement u is such that Aut = 0, i.e. 
Ut is the projection of the displacement u to the null space of the matrix A. The longitudinal 





Figure 1. Longitudinal and transverse waves in a simple quadratic lattice. 


component is a remaining orthogonal component of the displacement field and it is the 
projection of u to the row space of the matrix A. These projections have the following form [7] 

Ul = A?' Au, (6) 

Ut- = u — Ai^ {AA^)~^Au. (7) 

Indeed, Mt does not change the Voronoi cell volumes 

Aut- = Au — AAi^ {AA"^)~^ Au = 0 ( 8 ) 

and Ml is orthogonal to Mt 

m^Mt = AJ' {AAi^)~^ A{u — A^{AA^)~^Au) = 0. (9) 

3. The derivation of the matrix A 

Let us show how matrix A can be derived from geometry only. Let we shift only the atom j by 
the vector u (Figure]^. The vector u is small so the Voronoi cells after shifting have almost the 
same structure, but their facets are slightly shifted and rotated in space. The shifting u of the 
atom j can change the volume Vi of a nearest neighbor cell i. By definition the cells i and j are 
nearest neighbors if they have a common facet Sij. 

By definition of the Voronoi cell the facet lies in the plane Vij, which has the normal 
rvij = Tij/rij , where Yij = rj — r* is the vector connecting atoms i and j and go through the 
point bjj = (rj + rj)/2. Therefore the equation for this plane reads 

Uy • (r - bij) = 0. (10) 

After shifting of the atom j the new facet 5L lies in the plane "PL which satisfy the equation 

nij-(r-bV)=0 (11) 

where nV = I'ij/j'ij, + u, and bV = bj^ + u/2. The signed distance from an arbitrary 

point r to the plane PL is 

c^ii(r) = n7 • (b'j - r). (12) 

This distance has the sign if r lies on the same side to the plane PL as the atom i and the 
sign ” if r lies on the opposite side to the plane. If r lies in the plane PL then the distance 











































Figure 2. a) Two dimensional example of the Voronoi cells. The points B and C denote the 
position of hij and Cij respectively. In two dimensional case the point C is the center of the 
segment Sij. b) Dashed lines shows the Voronoi cells after shifting of the atom j by the vector 
u. The point B' denotes the position of h'-j. Green and red areas show the increasing and 
decreasing of the V respectively. 


d^j{r) is equal to 0 and we get the Eq. (Illl). In the linear approximation on u the change of the 
volume Vi is the integral of dij(r) over tEe surface of the initial facet Sij 


5Vi = 


T^ij ■ (b'j -r)dS = Sijtiii ■ (bE - c^), 
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where Sij and Cjj are correspondingly the area and the centroid of the facet Sij. In the linear 
approximation on u, the Eq. (13) reads 


XT/ 

OVi = - Y>ij ■ U, pij = Tj - Cij, 


(15) 




where we have taken into account that Cij ■ [hij — Cij) = 0 because Cjj G Sij. 

If all neighboring atoms shift, the change of the volume Vi can be written in the matrix form 


= X] 


,ja ^joL • 
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The nondiagonal elements of the matrix W we have already found above 

Wija = — (Pii)a, i + j- (17) 

fij 

The diagonal element Wija means the change of the ith Voronoi cell volume under shifting of 
the ith atom itself. Shifting of the all atoms by the same vector does not change the volumes of 
the Voronoi cells. Therefore 

V^ijci — ^ ^ kFjjQ, 

j¥=i 


( 18 ) 
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Figure 3. The decomposition of the total vibrational density of states to longitudinal 
and transverse components. The inset shows the relative number of the longitudinal modes 
gi^{ijj )/(red line). The relative number of the transverse modes grj.(Lo)/g{u}) = l — gi^{u})/g{uj) 
is shown by green hatching between red line and the value 1. 


After dividing by the volume we hnally get the matrix A 

It is notable that finite-elements methods also often use Voronoi cells and have a similar definition 
for the finite differences for the divergence operator [H]. 

4. Longitudinal and transverse DOS 

The total vibrational density of states can be decomposed to longitudinal and transverse 
components 


5(w) = - ^ j ) = 5l(w) + 5t(w), 
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where Uj is an eigenfrequency, u{ojj) is the corresponding eigenmode, Ui^{ojj) and are 

its transverse and longitudinal components defined by Eqs. and (0. The total vibrational 
density of states and its longitudinal and transverse components was calculated for the numerical 
model of a-Si consisting of = 32768 atoms by the Kernel polynomial method (KPM) [9]. The 
results are shown in the Figure 0 There are three clearly seen regions. The vibrations are 
mostly transverse below 7 THz and above 15 THz. In the frequency interval in between the 
vibrations have the longitudinal nature. 

The low number of longitudinal vibrations in the low-frequency region is explained by the 
Debye law, which gives 


£/l(w) ~ 


OJ 

■3 
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In our model the amorphous Si has the longitudinal sound velocity Cl = 7.96 km/s and the 
transverse sound velocity Ct = 3.85 km/s which gives the ratio /gT:{uj) = c|/2cl = 0.057 
for cj —7- 0. This value coincides with the inset in the Figure 

The domination of the longitudinal modes between 7 THz and 15 THz in amorphous silicon 
corresponds to the gap between the upper frequency of TA modes (7.5 THz) and the lower 
frequency of TO modes (13.9 THz) in crystalline Si [TO]. This frequency region in crystalline 
Si is totally occupied by LA and LO modes without a gap. In the same frequency region the 
vibrations of amorphous Si have a small transverse component (15-20%). Certainly, there are 
no optical phonons with well defined wavevector in amorphous silicon due to relatively strong 
disorder. However, the short-range order of vibrational modes in amorphous phase can be similar 
to that in crystal phase. At the same time, our definition of the longitudinal and transverse 
vibrations is local because it is based on the Voronoi cells, which depends on the neighbor atoms 
only. 

The Figure]^ shows a sharp change of the nature of vibrations at 7 THz from almost transverse 
to almost longitudinal. The longitudinal vibrations correspond to the stretching of the chemical 
bonds between atoms while the transverse vibrations correspond to the less rigid bond bending 
and rocking. Therefore longitudinal vibrations can transfer the vibrational energy much faster 
than transverse vibrations which leads to the sudden rise of the diffusivity at 7 THz observed 
in Refs. mm- Furthermore, our additional calculations show that the longitudinal modes still 
have a well-defined wavevector and a sound velocity up to ~ 13 THz unlike the transverse 
modes, which have the frequency of the loffe-Regel crossover wC/ ~ 4 THz. 

5. Conclusion 

We have shown that vibrations in amorphous solids can be decomposed to transverse and 
longitudinal components in a generalized sense. The transverse component does not change 
the volumes of Voronoi cells around atoms while the longitudinal component is the remaining 
orthogonal component changing this volume. We show that the rise of the diffusivity at 7 THz 
in amorphous silicon corresponds to the sharp change of the nature of vibrations from almost 
transverse to almost longitudinal ones having high sound velocity. This decomposition can be 
fruitful for the investigation of the correlation function and the diffusivity of transverse and 
longitudinal components of vibrations in amorphous solids independently from each other. 

The 3D-structure of several Voronoi cells in amorphous silicon is available online in 
supplementary materials. 
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